Blackbody radiation shift in a 43 Ca + ion optical frequency standard 
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Motivated by the prospect of an optical frequency standard based on 43 Ca + , we calculate the 
blackbody radiation (BBR) shift of the 4si/2 — 3ds/2 clock transition, which is a major component 
of the uncertainty budget. The calculations are based on the relativistic all-order single-double 
method where all single and double excitations of the Dirac-Fock wave function are included to all 
orders of perturbation theory. Additional calculations are conducted for the dominant contributions 
in order to evaluate some omitted high-order corrections and estimate the uncertainties of the final 
results. The BBR shift obtained for this transition is 0.38(1) Hz. The tensor polarizability of the 
3d 5 / 2 level is also calculated and its uncertainty is evaluated as well. Our results are compared with 
other calculations. 



PACS numbers: 06.30.Ft,32.10.Dk,31.25.-v,32.70.Cs 

The definition for the International System of Units 
(SI) of time, the second, is based on the microwave transi- 
tion between the two hyperfine levels of the ground state 
of 133 Cs. The advancements in experimental techniques 
such as laser frequency stabilization, atomic cooling and 
trapping, etc. have made the realization of the SI unit of 
time possible to 15 digits. The operation of atomic clocks 
is generally carried out at room temperature, whereas the 
definition of the second refers to the clock transition in 
an atom at absolute zero. This implies that the clock 
transition frequency should be corrected in practice for 
the effect of finite temperature of which the leading con- 
tributor is the blackbody radiation (BBR) shift. A num- 
ber of groups [l|, have accurately calculated the BBR 
shift at room temperature affecting the Cs microwave 
frequency standard. A significant further improvement 
in frequency standards is possible with the use of opti- 
cal transitions. The frequencies of feasible optical clock 
transitions are five orders of magnitude larger than the 
relevant microwave transition frequencies, thus making 
it theoretically possible to reach relative uncertainty of 
10 -18 . The ability to develop more precise optical fre- 
quency standards will open ways to improve Global Po- 
sitioning System (GPS) measurements and tracking of 
deep-space probes, perform more accurate measurements 
of the physical constants and tests of fundamental physics 
such as searches for nonlincarity of quantum mechanics, 
gravitational waves, etc. 

One of the promising schemes for the optical frequency 
standard with a single 43 Ca + ion trapped in a Paul trap 
was proposed by Champenois et al. [3[. A major advan- 
tage of using 43 Ca + ions is that the radiation required 
for cooling, repumping, and clock transition is easily pro- 
duced by non-bulky solid state or diode lasers. In this 
work, we calculate the shift in a 43 Ca + optical clock tran- 
sition due to blackbody radiation. The clock transition 
involved is 4si/ 2 — > 3d 5 / 2 F = 6 nip = 0. 

The accuracy of optical frequency standards is limited 



by the frequency shift in the clock transitions caused 
by the interaction of the ion with external fields. One 
such shift in the atomic levels is caused by external mag- 
netic fields (Zeeman shift). The first-order Zeeman shift 
can be eliminated by using isotopes with half-integer 
nuclear spin 0]. The second-order Zeeman shift de- 
pends on the choice of hyperfine sublevels involved in 
the optical transition. The particular transition scheme 
used to minimize second-order Zeeman shift in 43 Ca + is 
4si/2~^3<i5/2 F = 6 rap = 0. The effect of this shift on 
the clock transition is of the order of -0.09(9) Hz at 0.1 
fiT at room temperature [3]. There are some additional 
shifts in the clock transition frequency due to the second- 
order Doppler effect. These are on the order of mHz for 
ions at room temperature [H, H| . There have been a num- 
ber of accurate calculations of the Zeeman and Doppler 
shifts 3, 4]; however, there has been no accurate study of 
the BBR shift in this optical transition. As we show here, 
the BBR shift is the major component in the uncertainty 
budget of this frequency standard at room temperature. 

The frequency-dependent electric field E radiated by 
a black body at temperature T is given by the Planck's 
radiation law. The frequency shift of an ionic state due 
to such an electric field can be related to the states static 
scalar polarizability «o by 

Ai/ = -±(831.9 V/m) 2 (^.^ a (l + „), (1) 

where r/ is a small "dynamic" correction [5J] . We estimate 
that this dynamic correction is negligible with compari- 
son to the present 3 % accuracy of our calculations. The 
effect of the vector and tensor parts of the polarizability 
average out due to the isotropic nature of the electric field 
radiated by the black body. However, for completeness, 
we also present results for the static tensor polarizability 
of the 3d 5 / 2 state. 

The overall BBR shift of the clock transition frequency 
is calculated as the difference between the BBR shifts of 
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TABLE I: Contributions to the 4si/2 scalar (qo) static polar- 
izabilities in 43 Ca + and their uncertainties in units of a 3 ). The 
values of corresponding El matrix elements are given in eag. 



C out r ib ut ion 


(fcllDH4si /o) 

\ n 'l|- t -'ll'± £> l/2/ 




4Si/ 2 - 4pi /2 


2.898 


24.4(5) 


4S1/2 - 5pi/ 2 


0.076 


0.007 


4S1/2 - 6pi/ 2 


0.085 


0.007 


4S1/2 - 4p 3/2 


4.099 


48.4(1.0) 


4Si/2 - 5p 3/2 


0.091 


0.010 


4si/2 - 6p 3/2 


0.112 


0.012 


Q^corc 




3.25(17) 


Otail 




0.006(6) 


C*total 




76.1(1.1) 



the individual levels involved in the transition: 

A B br(4s -» 3d 5/2 ) = -- (ao(34/a) -«o(4s 1/2 )) 

x (831.9 V/m) 2 . (2) 

Therefore, the evaluation of the BBR shift involves ac- 
curate calculation of static scalar polarizabilities for the 
4si/ 2 ground and 3d 5 / 2 excited states. We use the rel- 
ativistic single-double (SD) all-order method to find the 
electric-dipole matrix elements used in the calculation of 
dominant contributions to the polarizability values. The 
relativistic all-order method is one of the most accurate 
methods currently being used in atomic structure calcu- 
lations. We refer the reader to Refs. [6|,l7j for the detailed 
description of this approach. 

The calculation of the static scalar polarizability of an 
atom can be separated into the calculations of the ionic 
core contribution a cor o and a valence contribution. The 
static ionic core polarizability values were calculated us- 
ing random-phase approximation (RPA) in Ref. [|| . The 
valence scalar and tensor polarizabilities of an atom in a 
state v are expressed as the sums over all intermediate 
states k allowed by the electric-dipole selection rules: 
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(k\\D\\v)' 
Ek — E v 



(3) 



5j v (2j v - 1) 



6(j v + l)(2j v + l)(2j v + 3) 



1/2 



(4) 



where (k \\D\\ v) are the reduced electric-dipole (El) ma- 
trix elements and Ei is the energy of a state i. In these 
equations, and hereafter, we use the conventional sys- 
tem of atomic units, a.u., in which e,m c , 47reo and the 
reduced Planck constant h have the numerical value 1. 
Polarizability in a.u. has the dimensions of volume, and 
its numerical values presented here are thus expressed in 
units of Qq, where ao ~ 0.052918 nm is the Bohr radius. 



TABLE II: Contributions to the 3d 5 / 2 scalar («o) and tensor 
(02) static polarizabilities in 43 Ca + and their uncertainties in 
units of a 3 ). The values of corresponding El matrix elements 
are given in eao. 



Contribution 


{k\\D\\3d 5/2 ) 


QO 


Ct2 


3rf 5 /2 


- 4p 3/2 


3.306 


22.78(25) 


-22.78(25) 


3^5/2 


- 5p 3/2 


0.148 


0.011(2) 


-0.011(2) 


oti 5/2 


— 6p 3 /2 


0.095 


0.004 


-0.004 


3^5/2 


~ 4/5/2 


0.516 


0.120(3) 


0.137(3) 


3^5/2 


- 5/5/2 


0.319 


0.039(2) 


0.044(2) 


3G?5/2 


- 6/5/2 


0.224 


0.018(1) 


0.020(1) 


3g?5/2 


~ 7/5/2 


0.169 


0.010 


0.011 


3^5/2 


- 8/5/2 


0.134 


0.006 


0.007 


3g?5/2 


- 9/5/2 


0.110 


0.004 


0.004 


3^5/2 - 


" IO/5/2 


0.093 


0.003 


0.003 


3(4/2 - 


- H/5/2 


0.079 


0.002 


0.002 


/2 




0.069 


0.002 


0.002 


3G?5/2 


- 4/7/2 


2.309 


2.392(53) 


-0.854(19) 


3c?5/2 


- 5/7/2 


1.428 


0.773(33) 


-0.276(12) 


3^5/2 


- 6/7/2 


1.000 


0.350(12) 


-0.125(4) 


3^5/2 


— 7 /7/2 


0.757 


0.191(7) 


-0.068(3) 


34/2 


— 8/7/2 


0.600 


0.117(4) 


-0.042(2) 


3G?5/2 


~ 9/7/2 


0.492 


0.077(3) 


-0.028(1) 


3c?5/2 - 


" 10/7/2 


0.414 


0.054(2) 


-0.019(1) 


3c?5/2 - 


" 11/7/2 


0.354 


0.039(1) 


-0.014(1) 


3g?5/2 - 


" 12/7/2 


0.308 


0.029(1) 


-0.011 


a 


core 




3.25(17) 








1.7(1.1) 


-0.5(3) 


Qftotal 




32.0(1.1) 


-24.5(4) 



The atomic units for a can be converted to SI units via 
a/h [Hz/(V/m) 2 ]=2.48832xl0~ 8 a [a.u.], where the con- 
version coefficient is 47reoag/ft, and Planck constant h is 
factored out. Experimental energies from Ref. @ have 
been used for the dominant contributions to the polariz- 
ability. 

We use the B-spline method to reduce infinite sums in 
Eqs. ^ and Q to a finite number of terms [l(| with 
70 splines of order 11 for each angular momentum con- 
strained to a spherical cavity with R = 220 a.u. The 
particular choice of such a large cavity ensures that all 
valence orbitals of interest are enclosed. In our case, the 
relevant excited states up to 12d states fit inside this cav- 
ity. We also separate the polarizability calculation into 
two parts, the main term ct ma in containing the first few 
contributions from the states that fit inside the cavity, 
with the remaining part designated by ataii- The main 
contribution is calculated using the all-order matrix el- 
ements obtained in the present work and experimental 
energies @. The tail contribution is evaluated in the 
Dirac-Fock (DF) approximation. We included as many 
states as practical in a ma in for the 3d 5 / 2 state in order to 
reduce the uncertainty in the remainder. 

Table Q] shows the contributions from the individual 
transitions to the ao of the 4s]/2 ground state of 43 Ca + . 
The main contributions are listed separately along with 
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the respective values of the electric-dipole matrix ele- 
ments. The tail contributions are grouped together as 
ataii- The total contribution of the As±/2 — 4pi/2 and 
4si/2 — 4p 3 / 2 transitions is overwhelmingly dominant. We 
use our ab initio SD all-order value for these matrix ele- 
ments based on the comparison of the all-order results for 
alkali-metal atoms with experiment [7j and assign these 
values 1 % uncertainty. 

The lifetimes of the Api/2 and 4f> 3 /2 states have been 
measured by Jin and Church [11] using a variant of 
the collinear laser-beam - ion-beam spectroscopy tech- 
nique in 1993. The experimental results, r(4p 1 / 2 ) = 
7.098(20) ns and T(4p 3/2 ) = 6.924(19) ns [n| were found 
to be in significant disagreement with perturbation the- 
ory calculations at that time. By 1993, similar discrep- 
ancies were found between other measurements carried 
out using the same technique and theoretical calcula- 
tions of light alkali-metal atom lifetimes. These dis- 
agreements in alkali-metal atoms were resolved a few 
years later when new measurements became available 
(see 0, EH and references therein). In all cases, new 
experimental results were found to be in excellent agree- 
ment with precise theoretical calculations. Recently, ex- 
cellent agreement was also found between experimen- 
tal and theoretical calculations of the 3d lifetimes in 
Ca + [l3j. Just as in the case of earlier comparison of 
theory and experiment for alkali- metal atoms, we find 
significant discrepancy between out present calculations 
yielding r(Api/ 2 ) = 6.875 ns, r(4p 3 / 2 ) = 6.686 ns and 
measurements of Ref. [ll|. It is interesting to make a 
direct comparison of the Apj — As theoretical and ex- 
perimental matrix elements. The contributions of the 
Ap — 3d transitions into the Ap lifetimes are small, as 
shown by comparison of our theory values for the Ein- 
stein A-coefficients A a f, for the transitions to the 4s 
state, A 4pi/2 4 S = 136.0 MHz and A ip3/2is = 139.7 MHz, 
and for the transitions to the 3dj states, A ipi/2 ?,d 3/2 = 
9.452 MHz, A 4p3/23d3/2 = 0.997 MHz, and A 4p3/2 3d 6/2 = 
8.877 MHz. Substituting the theoretical values for the 
A a i,, b — 3^3/2) 3^5/2, and experimental Apj lifetimes 
from [1 ll ] into the formula r a = l/X^^a&j we obtain 
the following values for the Apj — 4s reduced electric- 
dipole matrix elements: (4pi/ 2 ||d||4s} e:!::pt = 2.849(4) and 
(4p 1/2 \\d\\4s) expt = 4.023(6). While the uncertainties of 
these values include the uncertainties in the theoretical 
values of the Ap — 3d matrix elements, their contribu- 
tions are negligible. Our theoretical values differ from 
these experimental matrix elements by 1.7% and 1.9% 
for the Api/2 — As and 4p 3 / 2 —4s transitions, respectively. 
Our theoretical values for these transitions in K agree 
with experiment values within experimental uncertain- 
ties (0.13%). We note some differences between the con- 
tributions of the various correlation correction terms in 
the K and Ca + calculations of these matrix elements. 
Therefore, it would very interesting to see a new mea- 
surement of these matrix elements in Ca + to determine 
if the discrepancies between the theoretical and experi- 
mental values of the Ap lifetimes will be resolved as in 



TABLE III: Comparison of static scalar polarizabilities for 
4si/2 and 3d 5 / 2 states and blackbody radiation shift for the 
Asi/2 — 3^5/2 transition of 43 Ca + ion at T=300 K. The polar- 
izability values are in ajj and BBR shift is in Hz. 





Present work 


Ref. [3] 


Ref. [4] 


Ref. [14] 


ao(4s 1/2 ) 


76.1(1.1) 


76 


73 


70.89(15) 


ao(3d 5/2 ) 


32.0(1.1) 


31 


23 




BBR shift 


0.38(1) 


0.39(27) 


0.4 





the case of the alkali-metal atoms. 

Our all-order SD values are used for other main term 
contributions as well but their contribution is negligible. 
The uncertainty of the core contribution is taken to be 
5 % based on the comparison of the RPA approximation 
with experimental values for noble gases. 

Tabic HT1 shows the detailed breakdown of the contribu- 
tions to the 3d 5 / 2 state polarizabilities. In this case, the 
allowed transitions are to np^/2 and «/s/2.7/2 states. The 
two dominant contributions come from the 3c?5/2 — 4p 3 / 2 
and 3d 5 / 2 — 4/ 7 / 2 transitions. The sum of the contri- 
bution from the 3d 5 / 2 — nf transitions converges more 
slowly than the ns — np sums; therefore, we include in the 
main term the nine lowest 3d^/2~ n fj transitions for each 
j. In addition to the SD all-order calculations, we carry 
out semi-empirical scaling procedure for the two lowest 
3^5/2 — n P3/2 transitions and the three lowest 3d 5 / 2 — nf 
transitions to evaluate some classes of the correlation cor- 
rections omitted by the current all-order calculation. 

We also conduct additional more complicated ab ini- 
tio calculations that partially include triple excitation 
(SDpT) since those appear to be important for these 
types of transitions. The scaling procedure and SDpT 
calculations are described in Refs. 0, EU EH- We con- 
ducted a more detailed study of the uncertainties in the 
values of the final matrix elements for the first transi- 
tion in each sum, i.e. 3d$/2 — 4p 3 /2, id^/ 2 — Af 5 / 2 , and 
3^5/2 — 4/ 7 /2- They are determined as the maximum dif- 
ference between the final results and the ab initio and 
scaled SDpT values. The remaining uncertainties for the 
3d 5 /2~nf contributions are determined as the differences 
between the final values (either scaled SD or SDpT) and 
the ab initio SD values. The all-order calculation of such 
a large number of matrix elements and evaluation of their 
uncertainties makes the calculation of the 3c?5 /2 state po- 
larizability much more challenging than the calculation 
of the ground state polarizability. 

The tail contribution of the 3d 5 / 2 state is almost 5% of 
the total polarizability. Such a large tail leads to a large 
uncertainty in the total polarizability. In order to evalu- 
ate the uncertainty in the tail, we calculated the last few 
main terms using the DF approximation and compared 
the results with our all-order values. We found that the 
DF approximation overestimates the polarizability con- 
tributions by about 60 % relative to the all-order value. 
To improve our accuracy, we scaled the tail contribution 
with this ratio and took the difference of the DF and 
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scaled DF result to be the uncertainty of the tail value. 
The uncertainties from all terms are added in quadra- 
ture to obtain the uncertainties of the final polarizability 
values. 

We also list the contributions from various transitions 
to the 3g?5/2 tensor polarizability in Table HH The 3e?5/2 — 
4^3/2 transition also gives the dominant contribution to 
the tensor polarizability. The tail contribution is less 
important for a-i and the uncertainty in ai calculation is 
less than that of a$. 

We use the scalar polarizability values to evaluate the 
shift of a clock transition due to blackbody radiation at 
300 K. The final shift in the clock transition for a 43 Ca + 
ion is found to be (0.380 ± 0.013) Hz. The overall un- 
certainty of the result is dominated by the contribution 
of highly-excited nfj / 2 states to 3d 5 / 2 polarizability. 

Black body radiation also shifts the hyperfine split- 
tings of the ground and excited states P, although 
this is third-order effect and thus expected to be small. 
For completeness, we have calculated this effect which is 
proportional to the difference in the atomic polarizabil- 
ities between the two ground hyperfine structure states 
F of the atom. This is exactly the sort of calculation 
that is involved in the determination of the BBR shift in 
the microwave standard. The F-dependent polarizability 
can be written as the sum of three terms, each involv- 
ing a reduced hyperfine matrix element and two reduced 
electric-dipole matrix elements 0, [|| . We estimated the 
magnitude of the frequency shift of the ground state hy- 
perfine transition using DF values of these reduced ma- 
trix elements. For consistency, we used the same set of 
B-splines. Further details of the method can be found 
in Refs. [HQ- The frequency shift in the ground-state 
hyperfine transition of 43 Ca + at room temperature is cal- 



culated to be 4.2 x 10" 6 Hz. The magnitude of this third- 
order correction is approximately 10 -5 times lower than 
the second-order BBR shift and hence can be neglected. 
The effect of electric field on the frequency of the 3^/2 
state hyperfine transition is even smaller and; therefore, 
is not expected to affect the BBR shift of the clock tran- 
sition. 

In Table HTT1 we compare our polarizability and BBR 
shift results with previous theoretical studies. In the 
work done by Champenois et al. @, only first few tran- 
sitions were included in the Stark shift calculations. The 
required oscillator strengths were taken from the Kurucz 
database [lj] in Ref. Q and from Ref. [13 in Ref. [|. 
Our value for the ground state polarizability differs from 
that of Theodosiou [l4[ in part because of his omission of 
the core polarizability contribution in [l4j and his use of 
the Jin and Church [ll| values of the 4s — 4p matrix ele- 
ments. Our all-order values are in good agreement with 
those predicted by Champenois et al. [|[ and Kajita et 
al. 0]. We have significantly improved the accuracy of 
the BBR shift in comparison to these calculations by us- 
ing all-order values for the dominant terms and also by 
including the tail and core contributions. 

We have increased the accuracy of the polarizability 
values as compared to any other previous calculations by 
using a more accurate dipole-matrix element values cal- 
culated using relativistic all-order single-double method 
making the BBR shift value accurate enough to evaluate 
the performance of a Ca + single-ion frequency standard 
at room temperature. 

This research was performed under the sponsorship of 
the US Department of Commerce, National Institute of 
Standards and Technology. 
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